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Background: It is generally acknowledged that the time-dependent Hartree-Fock (TDHF) method provides a 
useful foundation for a fully microscopic many-body theory of low-energy heavy-ion reactions. The TDHF method 
is most widely known in nuclear physics in the small amplitude domain, where it provides a useful description of 
collective states, and is based on the mean-field formalism which has been a relatively successful approximation to 
the nuclear many-body problem. Currently, the TDHF theory is being widely used in the study of fusion excitation 
functions, fission, deep-inelastic scattering of heavy mass systems, while providing a natural foundation for many 
other studies. 

Purpose: With the advancement of computational power it is now possible to undertake TDHF calculations 
without any symmetry assumptions and incorporate the major strides made by the nuclear structure community 
in improving the energy density functionals used in these calculations. In particular, time-odd and tensor terms in 
these functionals are naturally present during the dynamical evolution, while being absent or minimally important 
for most static calculations. The parameters of these terms are determined by the requirement of Galilean 
invariance or local gauge invariance but their significance for the reaction dynamics have not been fully studied. 
This work addresses this question with emphasis on the tensor force. 

Method: The full version of the Skyrme force, including terms arising only from the Skyrme tensor force, is 
applied to the study of collisions within a completely symmetry-unrestricted TDHF implementation. 

Results: We examine the effect on fusion thresholds with and without the tensor force terms and find an effect 
on the fusion threshold energy of the order several MeV. Details of the distribution of the energy within terms in 
the energy density functional is also discussed. 

Conclusions: Terms in the energy density functional linked to the tensor force can play a non-negligible role in 
dynamic processes in nuclei, including those in which the terms do not affect ground state properties. 

PACS numbers: 21.60.Jz 


I. INTRODUCTION 


The extent to which the time-dependent mean field is 
a sufficiently good description of certain phenomena in 
nuclear dynamics is still an open question, as evinced 
partly by the ongoing explorations and developments in 
the area [1-14]. In particular, calculations of nuclear fu¬ 
sion in a time-dependent Hartree-Fock framework have 
developed over the years with increasingly relaxed sym¬ 
metry assumptions, and more sophisticated energy func¬ 
tionals in an attempt to settle this question. Approxi¬ 
mations of any type limit the number of degrees of free¬ 
dom accessible during a collision, and hence affect the 
nature and degree of dissipation. The understanding of 
the dissipative mechanisms in the TDHF theory is vital 
for establishing the region of validity of the mean-field ap¬ 
proximation and providing estimates for the importance 
of the mean-field effects at higher energies. In TDHF, the 
dissipation of the translational kinetic energy of the two 
ions is due to the collisions of single particle states with 
the walls of the time-dependent potential. This leads 
to the randomization of the motion characterized by the 
distribution of energy among all possible degrees of free¬ 
dom of the system. The complete equilibration of the 
translational kinetic energy among all possible degrees 
of freedom is commonly accepted as being the definition 


of fusion whereas the incomplete equilibration results in 
inelastic collisions. 

The inclusion of the spin-orbit interaction had a dra¬ 
matic effect on dissipation modes in heavy-ion colli¬ 
sions [15, 16]. The relaxation of any spatial symmetries 
gave rise to yet new modes [17]. The inclusion of all 
time-odd terms that arise from basic (not including ten¬ 
sor terms) Skyrme functionals was analyzed and found to 
have a noticeable effect on fusion properties [18] and in 
giant resonances [19]. These kinds of calculations help to 
pin down details of the nuclear energy density functional 
in ways complementary to studies of nuclear matter [20- 
23] and of the structure of finite nuclear structure [24-26]. 

In this work, we follow up previous studies by analyz¬ 
ing the effect of using Skyrme functionals which include 
the tensor part of the original Skyrme force, and we in¬ 
clude the most general terms in the density functional 
that arise, noting the interesting effects that have been 
seen in calculations of giant resonances when adding the 
tensor terms [27]. 

In Section H we briefly outline the Skyrme force and 
the resulting energy density functional used in our anal¬ 
ysis. Section HI presents TDHF calculations along with 
a discussion of the results, with a focus on the tensor 
terms of the functional. For comparison with previous 
work we concentrate on on collisions, where the 
tensor force has little effect on the static ground state 
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properties, so that in addition to comparison to previous 
studies, dynamic effects can be isolated. A concluding 
section (IV) follows the results. 


II. THE SKYRME FORCE AND THE ENERGY 
DENSITY FUNCTIONAL 


Skyrme’s interaction was originally posited as a zero- 
range low-momentum expansion of the effective inter¬ 
action in nuclear medium [28]. The original form for 
the two-body part of the potential was given as ti 2 = 
i5(ri — r 2 )t(k',k), where 


t(k', k) — to(l + xqP'^) + + k^) 

+ t 2 il + X 2 {P^ - h)k' ■ k 

5 

-f • k - -CTI • (T2fc^ -b conj.)) 

H— U{cri ■ k'(T2 ■ k -cti • (T2k' ■ k + conj.) 

2 3 

+ V(i{(Ti + (T2) ■ k'X k). ( 1 ) 


In the above equation, the third and fourth lines repre¬ 
sent the Skyrme tensor interaction, with T and U be¬ 
ing parameters to fit to data. In the first Hartree-Fock 
calculations using the Skyrme interaction, in 1972, the 
tensor terms were neglected [29]. This was reasonable, 
since only ground states of doubly-magic nuclei were cal¬ 
culated. The contribution of the tensor terms is, in that 
case, mainly to the spin-orbit splitting, which already has 
an adjustable parameter (V in (1)), and the data to which 
the force was htted did not demand extra parameters in 
the spin-orbit part. 

The effect of the tensor terms was studied later, in 
1977, again in doubly-magic nuclei [30], with mixed con¬ 
clusions as to their efficacy. Sporadically, parameter sets 
including tensor terms were explored, culminating in a 
recent resurgence in their study, motivated initially by 
the observed changing shell structure away from stability. 
The comprehensive paper by Lesinki et ah, gives a sum¬ 
mary of the history of the Skyrme tensor term, and we 
refer the reader there for a more complete account [24]. 

It is now routine to present the Skyrme interaction as 
an energy density functional (EDF). This is a more phys¬ 
ical approach given that the original three-body force has 
been generalized to a density-dependent two-body term 
with a fractional power of the density. Following Lesinski 
et al. [24], we present the full form of the Skyrme EDF, 
including terms which arise from the tensor force: 


£ = + cnpo]s^t + 
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where the following densities and currents are defined [24, 
31, 32] in terms of the density matrix in coordinate space 
for protons and neutrons (indicated by the subscript q): 

Pq{ra, r'a') = ]^Pq{r, -b ]^Sq{r, r') ■ (crVk)> (3) 

with 


Pqi'T, r') = ^ pq{ra, r'a') (4) 

<7 


and 


Sq{r, r') = ^ pq{ra, r'a'){a\cr\a). (5) 

< 7 < 7 ' 

The densities and currents needed are then defined in 
terms of the density pq{r,r') and spin density Sq{r,r') 
as 

Pqir) = Pqir,r')l^^, 

Sqir) = Sq{r,r')l^^, 

Tq{r) = V • Vpq{r,r')l^^, 

TqAr) = V-Vsqfr,r')l^^, (6) 

ff) = -^(V - V)pq{r,r') 

T—r' 

Jq,t,u{r)= -^(V^ - v;,)sg,^(r,r') 

r—r' 

FqAr)= i^(V^V', + V;,V.)s,..(r,r') 

- X ya - ya ^ 

where the Greek letter subscripts indicate Cartesian co¬ 
ordinates. From these densities one then defines the 
isoscalar (t = 0) and isovector (t = 1) densities and cur¬ 
rents found in (2) as 

Po{r) = pn{r)+pp{r) 

Pi{r) = p„{r) - pp{r), (7) 

and similarly of the other densities and currents. 
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In the version of the energy density functional pre¬ 
sented in (2), the contribution of the terms bilinear in 
the spin-current pseudotensor are presented in Carte¬ 
sian form. Some authors [32] use a pseudoscalar, vector 
and rank-2 pseudotensor representation instead. Denot¬ 


ing these as and respectively, the combi¬ 

nations of the pseudotensor that appear in the energy 
density functional can alternatively be expressed as [24]: 


1 

2 





t(0)" 


7-(2) t(2) 

(8) 
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- V 

(9) 
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We compute and present these terms in both ways both 
as a check of the implementation and to understand the 
ways in which the density functional distributes energy 
across different terms. 

We point out also, that in (2), the tensor force param¬ 
eters T and U wholly determine coupling constants 
and Cf, while contributing also, along with the central 
terms, to Cj and . Full details of the mapping be¬ 
tween Skyrme force parameterizations and energy density 
functional coefficients are well-documented [24]. We take 
all terms in (2), as defined by this mapping, except that 
we set = 0 since these terms are known to 

cause spin instabilities. The inclusion of the tensor force 
brings new terms in the density functional in to play, as 
well as modifying the strength of pre-existing terms. We 
note that a recent study analyzed the effect of the tensor 
parameters on the spin-orbit force alone (i.e. not includ¬ 
ing the “new” terms in the functional) in the context of 
heavy-ion collisions [33]. One does not need to make a 
necessary link between the coefficients in a force repre¬ 
sentation to the coefficients of the EDF. The latter can 
be considered as the free parameters to be fitted to data. 
It is not the purpose of the present paper to produce 
fits, however, but rather to evaluate the role of the terms 
that arise from the tensor force, and the effect on non¬ 
zero tensor force parameters on other terms in the EDF. 
We use existing Skyrme-l-tensor fits from the literature, 
motivated as follows. 

We choose SLy5 [35] as our test-bed, since there is a 
version available which has had the tensor forces added 
perturbatively to give improved single particle spectra in 
heavy isotopes without affecting the bulk properties [34]. 
It has also been extensively tested in TDHF collision cal¬ 
culations in which the full tensor was activated (as 
it was in the original SLy5 fit [35]) and time-odd terms 
studied [18]. Therefore, with this choice we will not have 
to re-fit a force, and are able to compare with the pre¬ 
vious complete calculations of collisions in TDHF, albeit 
with no tensor force. 

To test systematic properties of the tensor force pa¬ 
rameters, we also examine the T/J forces [24], in which 
the isovector and isoscalar parts of the tensor force are 


systematically varied, though fitted only through their 
contribution to the tensor. 


III. COLLISIONS 

Our calculations are performed on a Cartesian coor¬ 
dinate space grid with no symmetry assumptions, us¬ 
ing a version of the Sky3D code [37] with all time-odd 
and tensor terms included. In the first step, the ground 
state of is calculated with a damped gradient op¬ 
erator [38, 39]. Iterations continue until the variance in 
the single particle Hamiltonian is sufficiently small that 
the time-dependent calculation will be stable and the nu¬ 
cleus will translate without loss of energy on the grid [17]. 
We turn off any center of mass correction to be consistent 
between the static and dynamic calculations. Further de¬ 
tails of the set-up of our method can be found elsewhere 
[36] 


A. Upper Fusion Thresholds for head-on collisions 

The fusion window for heavy-ion collisions occurs when 
the nuclei have sufficient kinetic energy to overcome the 
Coulomb barrier, and sufficiently little that it can all be 
transferred to internal energy of the compound nucleus 
which then stays fused. Above the upper fusion thresh¬ 
old, the collision is deep-inelastic in nature and no fusion 
occurs. The fusion threshold can be rather sensitive to 
the energy density functional, particularly when different 
terms come in to play, which are not active in the ground 
states. We repeat previous calculations for SkM* [18], as 
a check, in which it is found that activating the terms 
(which are time-even) using the standard link between 
the Skyrme parameters and those in (2), a reduction in 
the fusion threshold of 6 MeV is found. We reproduced 
the size of the reduction, but found that activating the 
time-odd terms increased slightly the fusion threshold, 
bringing it slightly closer to the “basic” SkM* force, as 
shown in Table I 
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Force 

Threshold (MeV) 

SkM* (basic) 

77 

SkM* (inc. J^) 

71 

SkM* (full) 

73 

SLy5 (full) 

68 

SLy5t 

65 

ri2 

61 

ri4 

69 

T22 

64 

T24 

71 

T26 

82 

T42 

69 

TU 

79 

r46 

87 


TABLE I. Upper fusion threshold energies for the 
collision using various parameterizations of the Skyrme inter¬ 
action. 

One sees that the tensor terms, when added to an exist¬ 
ing parameterizations (i.e. SLySt compared to SLy5) can 
have a non-negligible effect of around 5% on the upper 
fusion threshold. The variation in results from the T/J 
parameterizations is more pronounced, yielding around 
a 25 MeV variation in the location of the upper fusion 
threshold. The lower threshold in each case is determined 
by the Coulomb barrier, so is equal between the forces. 
The fusion window, therefore, can be varied widely by 
the tensor part of the Skyrme interaction. 

B. EDF contributions in deep-inelastic collisions 

We examine the contribution from different parts of the 
Skyrme Energy Density Functional in the regime above 
the upper fusion threshold, so that the collision separates 
into parts before collision, during it, and afterwards when 
two nuclei re-emerge with internal excitation. For on 
at 100 MeV center of mass energy we look at SLy5 
and SLySt to examine the different contribution from the 
different terms. Figure 1 shows the contribution from the 

(split into those parts from the to and parameters), 
the s.T and the s.F terms. Respectively, these are terms 
which do not depend directly on the tensor coefficients 
(s^), are amended compared to the non-tensor case (s.T) 
and only come into play with non-zero tensor parameters 
{s.F). 

As shown in the two upper panels of Figure 1, the 
terms begin with zero contribution, as they should with 
two nuclei in their ground states, initialised with 
a Galilean velocity boost. Only when the nuclei begin 
to collide do these terms begin to differ from zero, as 
the collision process gives rise to regions of localised spin 
polarisation dynamically within the compound nucleus. 
Shortly after the collision, the contribution of the SLy5 
and SLy5t forces does not differ, as expected since the 



time [fm/c] time [fm/c] 



FIG. 1. (Color online) Contributions from terms involving 
some of the time-odd densities and currents to the total energy 
for 100 MeV collisions of on ^®0. 


coupling constants have not changed. Only later when 
the presence of the tensor terms has changed the over¬ 
all dynamics do the details of the contributions differ 
between the forces, owing to changes in the s density 
itself. 

A much more evident change in the dynamics is seen in 
the s.T and s.F terms. Without the tensor parameters, 
the s.F term is identically zero, while it is activated with 
a few hundred keV of the available energy during collision 
via coupling to the tensor parameters. The s.T terms, on 
the other hand, have couplings combining the tensor and 
surface terms of the Skyrme interaction that conspire to 
much reduce the role of this term when the tensor terms 
are activated in SLySt, compared to SLy5. 

Both the terms s.T and s.F are linked by Galilean 
invariance to terms arising from bilinear couplings of the 
spin-current tensor J, as given by those terms in (2) 
which share the same parameter. 

As a check of our approach, and also to understand the 
underlying dynamics and interplay between the terms in 
the functional, we evaluate the contributions due to the 
terms both in their Cartesian form as given in (2), and 
in their coupled form given in (8) and (9). Figure 2 shows 
these contributions. The column on the left shows the en¬ 
ergy contribution from the (pseudo)scalar-, vector- and 
(pseudo)tensor-decomposed form of while the right 
hand side shows those as they appear in (2). The first 
three plots in the right column are the diagonal, antisym¬ 
metric and symmetric terms respectively. As expected, 
the first and third terms are identically zero in the case 
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FIG. 2. (Color online) Contributions to the total energy from 
the terms following the Cartesian decomposition as in the 
paper by Lesinski [24] (right column) or following the coupled 
form given by Perlinska [32] (left column). The bottom frames 
show the total contribution in both cases. The case is 100 
MeV collisions of on ^®0. 



FIG. 3. (Color online) Contribution to the total energy from 
the terms for a 6 = 6.65fm collision, which just fuses, at 
34 MeV as a function of force. 


nucleus. In Figure 3 we plot the contribution to the total 
energy of the system as described, at 6 = 6.65 fm, from 
the total terms as a function of the force parameters. 
It is seen that the perturbative addition of the tensor 
terms to the SLy5t forces results in a different sign of 
the contribution to the mean field. The results from 
the TIJ forces span a range of contribution, with up to 
several MeV of energy being stored in this term at times 
during the collision process. The terms are not identi¬ 
cally zero at t = 0, but they are greatly excited during the 
collision process. It is clear that the different behaviour 
of different tensor parameterisations can appear in such 
a dynamic situation, while being much less evident in the 
ground state (to which the force is fitted). 

We mention in this section that the impact parameter 
dividing those configurations which fuse from those which 
don’t was rather insensitive to the tensor parameter set, 
at least for this center of mass energy of 34 MeV. There is 
thus little effect on the cross-section. The study of Dai et 
al. [33] in which the tensor contribution to the spin-orbit 
interaction was evaluated, a variation in the cross-section 
is seen at 70.5 MeV. This is consistent with the fact that 
we see a resonable spread of upper thresholds around this 
energy. 


of SLy5 since they appear only multiplied by the tensor 
coefficients. The bottom frames of figure 2 show the total 
contribution when calculated by both approaches. They 
are identical, as they must be. 


C. Off-axis collisions 

For a series of different parameterisations, we per¬ 
formed calculations of on at a center-of-mass 
energy of 34 MeV and impact parameter close to a graz¬ 
ing impact, to validate against previous work [18] for the 
tensor-less forces and to understand the dynamics in the 


IV. CONCLUSIONS 

We have performed time-dependent energy density 
functional calculations of heavy ion collisions using 
^®0-|-^®0 as a test case. We have included in the density 
functional all terms that arise when deriving the func¬ 
tional from a Skyrme force including tensor terms, and 
used existing Skyrme-l-tensor parameterisations to assess 
the effect of the tensor terms within this framework. It 
is found that the size of the fusion window can vary as 
the tensor force varies, owing to a movement in the up¬ 
per fusion threshold. Contributions from different terms 
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in the functionals were analysed. Large variations in the 
dominant terms, which could be either attractive or 
repulsive, were found for different tensor-dependent func¬ 
tionals. 


Terms in the density functional which arise from the 
tensor term, which play a minor role in fitting ground 
state data, can have a significant effect in dynamical 
properties. Thus there is scope for adjustment of func¬ 
tionals (or forces) in light of new data on dynamical pro¬ 
cesses while retaining good fits to ground state data. 
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